Foundations Series • From Euler to Hawking

Gauss & Maxwell: The Foundations of Quantum Mechanics and String Theory

Carl Friedrich Gauss gave us the mathematics of curvature and the Gaussian integral. James Clerk Maxwell unified electricity and magnetism into field theory. Together, their ideas are the scaffolding for Feynman's path integrals, Hilbert's quantum states, and the vibrating worldsheets of string theory.

Gaussian Integrals Gauge Theory Path Integrals Gauss-Bonnet
Carl Friedrich Gauss portrait
1777–1855
James Clerk Maxwell portrait
1831–1879
Maxwell in Vacuum
$$ \begin{aligned} \nabla \cdot \mathbf{E} &= 0 \qquad \nabla \cdot \mathbf{B} = 0 \\ \nabla \times \mathbf{E} &= -\frac{\partial \mathbf{B}}{\partial t} \\ \nabla \times \mathbf{B} &= \mu_0\varepsilon_0\frac{\partial \mathbf{E}}{\partial t} \end{aligned} $$
1

Gauss's work is everywhere in modern physics. The Gaussian integral is the only integral we can do exactly in infinite dimensions—making quantum field theory possible.

$$ \int_{-\infty}^{\infty} e^{-a x^2}\,dx = \sqrt{\frac{\pi}{a}}, \quad \Re(a)>0 $$

The normalized Gaussian defines probability and quantum uncertainty:

$$ \frac{1}{\sqrt{2\pi\sigma^2}}e^{-(x-\mu)^2/2\sigma^2}, \quad \int_{-\infty}^{\infty} |\psi|^2 dx = 1 $$

Gauss's Law for electricity connects flux to charge—later a constraint in quantum gauge theory:

$$ \oint_{\partial V} \mathbf{E}\cdot d\mathbf{A} = \frac{Q_{\text{enc}}}{\varepsilon_0} $$
$$ \nabla\cdot\mathbf{E} = \frac{\rho}{\varepsilon_0} $$

For surfaces, Gaussian curvature $K$ and the Gauss-Bonnet theorem link local geometry to global topology—central to string worldsheets:

$$ \int_M K\,dA = 2\pi\chi(M), \qquad \chi = 2-2g $$

Vary $a$ in $e^{-ax^2}$. The area is always $\sqrt{\pi/a}$—the cornerstone of all Gaussian path integrals.

Integral = 1.772

This exact result lets Feynman and Hawking evaluate $\int \mathcal{D}x\, e^{iS/\hbar}$ for quadratic actions.

2

Maxwell unified electricity and magnetism. His equations imply light is an electromagnetic wave with $c = 1/\sqrt{\mu_0\varepsilon_0}$.

Modern form uses the gauge potential $A_\mu = (\phi/c, \mathbf{A})$ and field tensor:

$$ F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu, \quad \partial_\mu F^{\mu\nu}= \mu_0 J^\nu $$

Gauge freedom $A_\mu \to A_\mu + \partial_\mu \Lambda$ leaves $F_{\mu\nu}$ invariant. This $U(1)$ symmetry is the template for the Standard Model and for string theory's Chan-Paton factors.

Following Euler’s variational methods and Hilbert’s operator formalism, Maxwell’s theory becomes the first quantum field theory upon quantization.

Flux through any sphere around a charge is $Q/\varepsilon_0$, independent of radius.

Q = 1.0 e
r = 1.0
∮E·dA = 1.13e11 V·m
3

Quantum states live in Hilbert space. The Gaussian wave packet saturates Heisenberg:

$$ \psi(x,0) = \left(\frac{2a}{\pi}\right)^{1/4} e^{-a(x-x_0)^2} e^{ip_0 x/\hbar}, \quad \Delta x\,\Delta p = \frac{\hbar}{2} $$

Feynman's path integral is defined by Gaussian integrals. For a free particle:

$$ K(x_f,t;x_i,0)=\int \mathcal{D}x\, e^{iS/\hbar} = \sqrt{\frac{m}{2\pi i\hbar t}}\,e^{im(x_f-x_i)^2/2\hbar t} $$

All evaluations reduce to $\int e^{-a x^2 + bx}dx$. Gauss makes quantum mechanics computable.

Gauss's law as constraint: Physical states satisfy

$$ (\nabla\cdot\hat{\mathbf{E}} - \hat{\rho})|\psi_{\text{phys}}\rangle = 0 $$

removing unphysical polarizations—exactly as Dirac constrained Hilbert spaces.

Wider initial packets spread slower. This is quantum diffusion.

σ₀ = 0.50
Δx(t) = 0.50, Δp = 1.00
4

String theory’s Polyakov path integral is Gaussian in the embedding fields $X^\mu$:

$$ Z = \int \mathcal{D}X\,\mathcal{D}h\, e^{-S_P},\quad S_P = \frac{1}{4\pi\alpha'}\int d^2\sigma\sqrt{h}\,h^{ab}\partial_a X^\mu \partial_b X_\mu $$

Evaluating it gives the propagator $\langle X(z)X(w)\rangle \sim -\alpha'\ln|z-w|$, again a Gaussian integral.

The action contains the Gauss-Bonnet term:

$$ S \supset \lambda \int d^2\sigma\sqrt{h}\,R = 4\pi\lambda\chi = 4\pi\lambda(2-2g) $$

Thus amplitudes sum over genus $g$ with weight $g_s^{2g-2}$—topology controls quantum corrections.

Maxwell from strings: Open string massless mode is $A_\mu$. Its low-energy action is Born-Infeld, reducing to

$$ S_{\text{DBI}} \to -\frac{1}{4}\int F_{\mu\nu}F^{\mu\nu} d^4x $$

D-branes source RR fields, obeying Gauss's law: $\oint_{S^{8-p}} *F_{p+2}=Q_{\text{Dp}}$.

Gauss-Bonnet: $\chi = 2-2g$ determines string coupling weight.

χ = 2
∫K dA =
weight ∼ gₛ⁻²
5

Gauss → Path Integrals

The Gaussian integral → finite-dimensional Gaussian integrals → infinite-dimensional path integrals. Hawking's Euclidean quantum gravity uses the same $\int e^{-S_E}$.

Maxwell → Gauge Theory

$U(1)$ gauge symmetry generalizes to $SU(3)\times SU(2)\times U(1)$ in the Standard Model, and to non-abelian gauge fields on D-branes in string theory.

Gauss-Bonnet → Topology

Gauss's curvature theorem controls string loop expansion. Euler characteristic, refined by Riemann and Chern, classifies Calabi-Yau compactifications.

Building on Euler’s $e^{i\pi}+1=0$ and Hilbert’s operator algebras, Gauss and Maxwell provide the computational engine and physical principle that Hawking later used for black hole radiation.
6

Run these in Octave to verify the mathematics Gauss and Maxwell built. All scripts are self-contained and use only core functions.

Compare analytic $\sqrt{\pi/a}$ with numerical quadrature.

% Gaussian integral
a = 2;
I_analytic = sqrt(pi/a);
I_numeric = integral(@(x) exp(-a*x.^2), -Inf, Inf);
fprintf('Analytic: %.10f, Numeric: %.10f, Error: %.2e\n', ...
        I_analytic, I_numeric, abs(I_analytic-I_numeric));

x = linspace(-3,3,400);
plot(x, exp(-a*x.^2), 'LineWidth',2); grid on;
title('Gaussian e^{-ax^2}'); xlabel('x');

Initial minimum-uncertainty packet $|\psi|^2$.

% Wave packet at t=0
x = linspace(-5,5,500);
sigma0 = 0.5;
psi = (2*pi*sigma0^2)^(-0.25) * exp(-x.^2/(4*sigma0^2));
prob = abs(psi).^2;

plot(x, prob, 'LineWidth',2); grid on;
xlabel('x'); ylabel('|\psi|^2');
title(sprintf('Gaussian packet, \\sigma_0 = %.2f', sigma0));
% Uncertainty: Delta_x * Delta_p = hbar/2

Flux through sphere independent of radius.

% Gauss's law
Q = 1.602e-19;          % elementary charge
eps0 = 8.854187817e-12;
r = logspace(-9,-7,5);   % various radii
E = Q ./ (4*pi*eps0*r.^2);
flux = 4*pi*r.^2 .* E;  % = Q/eps0

fprintf('Q/eps0 = %.3e\n', Q/eps0);
disp([r' flux']);  % constant column

Plane wave solution $E = E_0\sin(kx-\omega t)$.

% EM wave
c = 299792458; E0 = 1;
k = 2*pi; omega = c*k;
x = linspace(0,2,500);
t = 0;
E = E0*sin(k*x - omega*t);
B = E0/c * sin(k*x - omega*t);

plot(x,E,'r',x,B*3e8,'b','LineWidth',1.5); grid on;
legend('E (V/m)','cB (V/m)'); xlabel('x (m)');
title('Maxwell wave satisfies \nabla^2E = \mu_0\epsilon_0 \partial_t^2E');

Verify $\int K dA = 4\pi$ for $S^2$.

% Gauss-Bonnet sphere
R = 1;
K = 1/R^2;                % Gaussian curvature
area = 4*pi*R^2;
integral_K = K * area;
chi = integral_K/(2*pi);  % Euler characteristic

fprintf('∫ K dA = %.5f = 2πχ, χ = %.1f\n', integral_K, chi);
% For genus g: chi = 2-2g, torus g=1 gives 0